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A  thermoelectric  element  made  of  p-  and  n-type  semiconductor  plates  bonded  onto  a  highly  thermal  and 
electrical  conducting  inter-connector  material  with  an  integrated  flow  channel  can  be  treated  as  an  inte¬ 
grated  thermoelectric  device  (iTED).  The  performance  of  an  iTED  with  multiple  elements  connected  elec¬ 
trically  in  series  and  thermally  in  parallel  has  been  investigated  using  numerical  simulations.  The  top  and 
bottom  surfaces  of  the  device  are  subjected  to  a  constant  cold  temperature  while  the  inter-connector 
channel  walls  are  exposed  to  a  hot  fluid.  The  thermoelectric-hydraulic  behavior  of  an  iTED  is  analyzed 
in  terms  of  heat  input,  power  output,  conversion  efficiency,  produced  electric  current,  Ohmic  and  Seebeck 
voltages,  and  pressure  drops  for  various  hot  fluid  flow  rates  Re  and  inlet  temperatures  Tin,  thermoelectric 
material  sizes  d,  and  number  of  modules  N.  For  a  single  module  iTED  with  fixed  d  and  Tin  values,  the 
power  output  and  efficiency  are  increased  five-  and  twofold,  respectively  at  Re  =  500  when  compared 
with  the  values  of  Re  =  100.  For  given  Re  and  d  values,  increasing  Tin  resulted  in  enhanced  device  perfor¬ 
mance.  Furthermore,  increasing  d  increased  internal  resistance  and  resulted  in  a  decrease  of  heat  input. 
The  influence  of  d  on  power  output  is  phenomenal;  for  a  given  set  of  geometric  and  thermal  boundary 
conditions,  there  exists  an  optimum  d  where  a  maximum  power  output  is  achieved.  The  addition  of  mod¬ 
ules  N  resulted  in  a  significant  improvement  in  power  output  and  a  reduction  in  produced  electric  current 
and  efficiency.  For  instance,  device  with  N  =  5  showed  more  than  a  twofold  increase  in  power  output  and 
nearly  a  33%  reduction  in  both  efficiency  and  electric  current  when  compared  to  N  =  1  values. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

There  is  an  ever-increasing  amount  of  green-house  gas  and 
waste-heat  released  into  the  atmosphere  from  the  fossil-fuel 
power  generation  plants,  automobiles,  and  industrial  heating  or 
cooling  systems  in  response  to  continual  energy  demands.  Approx¬ 
imately  two  -thirds  of  the  supplied  energy  into  these  systems  is 
rejected  as  a  waste-heat  to  the  surroundings.  There  is  an  urgent 
need  to  explore  novel,  environmentally-friendly  technologies  that 
can  replace  or  improve  the  performance  of  the  existing  systems. 
Solid-state  thermoelectric  devices  are  a  viable  technology  for 
recovering  waste-heat  and  convert  it  into  electricity  while  mitigat¬ 
ing  the  emission  of  green-house  gases. 

Thermoelectric  devices  (TEDs)  are  constructed  by  joining  two 
different  electrically  and  thermally  conductive  materials  at  a 
junction.  Using  the  Seebeck  effect,  TEDs  work  as  electric  power 
generators  when  the  two  material  junctions  are  exposed  to  a 
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temperature  differential.  Similarly,  TEDs  act  as  a  refrigerators  via 
the  Peltier  effect  when  an  electric  current  is  applied  across  the  ter¬ 
minals,  creating  a  temperature  differential  at  the  material  junction 
[1].  However,  the  current  thermoelectric  materials  with  a  figure  of 
merit  «1.5  achieve  thermal  conversion  efficiencies  of  5-15%  and 
coefficients  of  performance  (COP)  of  0.5-1.  Due  to  their  scalable, 
reliable,  stable,  compact  and  noise  free  operation,  TEDs  are  suitable 
in  novel  applications  such  as  waste-heat  recovery  from  exhaust 
streams  and  other  low-grade  heat  sources,  electric  power  genera¬ 
tion  for  remote  radio  and  satellite  stations,  pocket  electronics, 
bio-thermal  batteries  to  power  pacemakers,  localized  cooling  in 
electronic  components  and  space  cooling  in  automobile  seats. 

The  efficiency  of  TEDs  has  been  increased  via  the  methods  of 
nano-structuring  and  fabrication  [2-4],  novel  designs  [5-10]  and 
use  of  new  bulk  materials  [1].  Caillat  et  al.  [5]  developed  a  seg¬ 
mented  TED  using  novel  p-  and  n-type  materials,  and  achieved  a 
conversion  efficiency  of  15%.  El-Genk  et  al.  [6  reported  peak  effi¬ 
ciencies  of  16.69%  and  7.4%  respectively  for  skutterudite  and  SiGe 
segmented  TEDs.  Further,  Punnachaiya  et  al.  [7]  studied  cascaded 
TEDs  and  showed  a  low  conversion  efficiency  of  0.47%  with 
Th  =  96  °C  and  temperature  differential  ( Th-Tc )  of  25  °C.  Liang 
et  al.  [8]  investigated  the  performance  of  a  multistage  TED  and 
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Nomenclature 


A  cross-sectional  area,  mm2 

cp  specific  heat  of  fluid,  Jkg-1 K_1 

d  size  of  semiconductor  material,  mm 

D  depth  of  the  thermoelectric  leg,  mm 

Dh  hydraulic  diameter  of  main  flow  channel  ^H-ld+D) 
mm 

H  height  of  the  leg,  mm 

/  electric  current,  A 

J  electric  current  density  A  m-2 

J  current  density  vector 

k  thermal  conductivity,  Wm'1  K1 

L  distance  between  the  legs,  mm 

N  number  of  thermoelectric  modules 

P  pressure,  N  m-2 

P0  power  output  (I2RL),  W 

Q.  heat  transfer,  W 

R  electrical  resistance,  Q 

Re  Reynolds  number  ( pUDhlp ) 

t  connector  thickness,  mm 

T  temperature,  K 

U  inlet  velocity,  m  s-1 

v  velocity  vector 

V  voltage,  V 

u ,  v,  w  velocities  in  x,  y,  z  directions,  m  s-1 


W  width,  mm 

x,  y,  z  coordinates,  mm 

Greek  symbols 

a  Seebeck  coefficient  or  thermopower,  V  K_1 

p  dynamic  viscosity,  N  s  m-2 

p  electrical  resistivity,  Q  m 

pf  density  of  fluid,  kg  m-3 

rj  thermoelectric  conversion  efficiency,  dimensionless 

Subscripts 

c  cold  wall/conductor 

in  inlet 

/  fluid 

h  hot  wall 

i  internal/integrated 

ic  inter-connector 

L  load 

n  n-type  semiconductor 

0  Ohmic  potential 

p  p- type  semiconductor 

S  Seebeck  potential 

£  direction  normal  to  the  surface 


they  demonstrated  that  the  thermal  contact  resistance  between 
the  TED  module  and  heat  source  or  sink  plays  a  substantial  role 
in  the  power  output.  Recently,  Crane  et  al.  9]  built  a  full-scale 
cylindrical-shaped  thermoelectric  generator  using  segmented  and 
high-power  density  elements  and  produced  a  power  output  of 
608  W. 

Using  analytical  solutions  and  experiments,  Gou  et  al.  [11] 
showed  that  expanding  the  heat  sink  surface  area  and  enhancing 
the  cold  side  heat  transfer  in  proper  ranges  can  have  significant  ef¬ 
fects  on  TEDs  performance  when  compared  to  increasing  the 
waste-heat  temperature  and  addition  of  modules  in  series.  Further, 
Hasio  et  al.  [12]  investigated  TED  applied  to  waste-heat  recovery 
from  an  automobile  engine  and  showed  that  the  TED  performs  bet¬ 
ter  on  the  exhaust  pipe  than  radiator.  Even  though  electric  and 
temperature  fields  in  TEDs  are  multidimensional,  most  of  the 
works  in  literature  were  conducted  with  one-dimensional  analyti¬ 
cal  solutions  for  simplified  cases.  However,  for  accurate  system 
design  and  optimization  of  TEDs,  the  three-dimensional  (3D) 
numerical  studies  have  been  conducted  in  the  articles  [13- 
20,22].  Further,  Karri  et  al.  [21]  investigated  the  performance  of 
various  modeling  approaches  for  thermoelectric  elements  and 
their  merits  and  demerits  under  given  geometric  and  thermal 
boundary  conditions. 

Harris  et  al.  [13]  studied  the  influence  of  inert  gas  and  insulat¬ 
ing  materials,  and  interface  contact  resistances  on  TED  perfor¬ 
mance  using  finite  volume  numerical  methods.  Considering  the 
convection  and  radiation  effects,  Ziolkowski  et  al.  [14]  studied 
the  performance  of  TED  in  ANSYS  software  for  various  pellet  aspect 
ratios  and  contact  resistances.  Gould  et  al.  [22]  simulated  the  TED 
configured  for  low  power  generation  using  TCAD  package.  In  the 
articles  [15,16],  the  authors  used  numerical  simulations  to  study 
the  optimum  geometries  of  TEDs  to  achieve  a  maximum  conver¬ 
sion  efficiency.  Furthermore,  the  researchers  [17,19]  proposed 
and  implemented  a  3D  numerical  model  for  thermoelectric  gener¬ 
ators  in  FLUENT  UDS  environment  and  their  model  accounts  for  all 
temperature-dependent  properties  of  materials  and  non-linear 
fluid-thermal-electric  multi-physics  coupled  effects.  Recently,  by 


coupling  both  temperature  and  electric  potential  fields  Wang 
et  al.  [18,20]  studied  the  steady  and  transient  response  of  TEDs 
using  3D  multi-physics  models. 

It  is  observed  from  literature  that  the  conventional  TED  designs 
applied  to  waste-heat  recovery  induce  large  thermal  resistance  be¬ 
tween  the  working  fluid  and  the  thermoelectric  junctions  via  the 
heat  exchanger,  ceramic  plate,  and  the  interface  materials,  and  also 
they  require  great  amounts  of  semiconductor  materials  for  module 
construction  due  to  the  demand  of  constant  leg  heights  for  fabrica¬ 
tion  ease.  Keeping  this  in  mind,  the  authors  proposed  an  integrated 
thermoelectric  device  (iTED)  [23,24]  where  the  hot-side  inter-con¬ 
nector  is  re-designed  to  incorporate  an  integrated  heat  exchanger. 
By  doing  so,  this  novel  design  reduces  the  thermal  resistances 
attributed  to  the  heat  sinks  and  ceramic  plates,  which  are  present 
in  convectional  design.  Further,  by  adjusting  the  height  of  the  in¬ 
ter-connector  heat  exchanger,  this  design  would  help  in  using  dif¬ 
ferent  semiconductor  element  heights  for  achieving  maximum 
power  output  and  efficiency  while  keeping  the  height  of  the  mod¬ 
ule  invariant.  In  essence,  this  design  enhances  reliability  and  per¬ 
formance  of  TED.  In  this  study,  using  numerical  methods  the 
thermoelectric-hydraulic  performance  of  such  a  iTED  with  multi¬ 
ple  modules  connected  electrically  in  series  and  thermally  parallel 
is  investigated.  The  effects  of  heat  exchanger  hot  fluid  flow  rates 
and  inlet  temperatures,  thermoelectric  element  sizes,  and  the 
number  of  modules  on  the  performance  of  an  iTED  applied  to 
waste-heat  recovery  under  steady  state  conditions  (the  perfor¬ 
mance  variation  with  respect  to  the  time  is  zero)  are  studied  in 
detail. 


2.  Geometry,  governing  equations,  and  boundary  conditions 

The  schematic  of  the  three-dimensional  single-stage  integrated 
thermoelectric  device  (iTED)  being  investigated  is  shown  in  Fig.  la. 
The  device  consists  of  two  vertical  legs  connected  electrically  in 
series  and  thermally  in  parallel  via  connectors  made  of  highly  elec¬ 
trical  conducting  material.  Each  thermoelectric  leg  is  constructed 
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(b) 


Fig.  1.  Schematic  of  (a)  single-stage  integrated  thermoelectric  device  and  (b)  mesh. 


by  inserting  a  copper  inter-connector  with  an  integrated  flow 
channel  acting  as  a  heat  exchanger  between  n-  and  p- type  semi¬ 
conductor  slices.  A  square  cross-section  W  x  D  and  height  H  is  con¬ 
sidered  for  each  leg  and  the  legs  are  separated  by  a  distance  L  as 
shown  in  Fig.  la.  The  inter-connector  channel  of  rectangular  shape 
with  size  0.6D  x  0.8(L  -2 d)  x  W  is  taken  for  fluid  flow.  A  hot  fluid 
with  constant  temperature  Tin  and  uniform  velocity  U  enters  the 
main  flow  channel  of  invariant  cross-sectional  area  D  x  (I  -  2d) 
and  passes  through  inter-connector  channel  and  leaves  at  the  exit 
of  the  fluid  domain  (Fig.  la).  The  channel  walls  and  the  up  and 
down  stream  side  surfaces  of  inter-connector  are  solely  responsi¬ 
ble  for  the  heat  transfer  from  the  hot  fluid  to  the  thermoelectric 
elements.  The  top  surface  of  the  upper  connector  and  the  lower 
surface  of  the  bottom  connectors  are  maintained  at  a  constant  cold 
temperature  Tc.  The  remaining  surfaces  of  the  iTED  and  the  main 
flow  channel  walls  are  kept  adiabatic. 

The  iTED  generates  electrical  power  through  a  temperature  dif¬ 
ferential  ( Th  -  Tc )  present  between  the  inter-connector  channel 
walls  and  the  cold  surface  via  the  Seebeck  effect  when  it  is  con¬ 
nected  to  a  load  resistance  RL.  The  fluid  flow  and  the  heat  transport 
in  the  fluid  domain  are  governed  by  the  conservation  of  mass, 
momentum,  and  energy  equations  while  subjected  to  the  following 
assumptions:  the  flow  is  steady,  laminar  and  incompressible,  and 
the  thermophysical  properties  of  fluid  are  kept  constant.  Thermo¬ 
electric  materials  are  heterogeneous  and  isotropic.  Furthermore, 
the  current  flow  and  the  heat  transport  in  thermoelectric  materials 
are  governed  by  the  continuity  of  current  density  and  the  energy 
equations,  respectively. 

The  set  of  partial  differential  equations  governing  the 
incompressible  fluid  flow  and  heat  transport  in  the  fluid  do¬ 
main  and  the  thermoelectric  phenomenon  in  the  conductor 
and  semiconductor  materials  under  steady-state  conditions 
are  written  as: 


•  Continuity  equations: 

-  Fluid  region: 

V.v  =  0.  (1) 

-  Conductor  and  semiconductor  material  [1]: 

V.J  =  0.  (2) 

•  Momentum  equation: 

Pf(v.Vv)  =  VP  +  y(iV2v.  (3) 

•  Energy  equations: 

-  Fluid  region: 

(PCp)f{v-S7T)  =  V.(kfVT).  (4) 

-  Conductor  material: 

V.(kVT)  +  pj2  =  0.  (5) 

-  Semiconductor  material  [1]: 

V.(kVT)  +  pi2  -  TJ.  [(Va)T  +  (J0  VT]  =  0.  (6) 


The  second  and  third  terms  in  the  left  hand  side  of  Eq.  (6)  rep¬ 
resent  the  Joule  heating  and  Peltier  and  Thomson  effects,  respec¬ 
tively.  The  Thomson  effect  as  compared  to  Joule  heating  and 
Peltier  effect  does  not  contribute  significantly  to  the  performance 
of  thermoelectric  material.  However,  it  has  been  included  in  the 
analysis  for  the  completeness. 

Using  the  non-Ohmic  current-voltage  [25]  relation,  the  total 
electric  potential  is  calculated  as 

W  =  Wo  +  Ws  =  -pi  -  aVT 


(7) 
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In  Eq.  (7),  the  first  term  on  the  right  hand  side  is  the  electrostatic 
potential  distribution  due  to  current  flowing  in  the  material  and 
the  second  term  comprises  the  Seebeck  electric  potential  distribu¬ 
tion  via  the  temperature  differential  created  in  the  thermoelectric 
material. 

With  respect  to  the  geometry  shown  in  Fig.  la,  the  associated 
thermal  and  electrical  boundary  conditions  for  solving 
Eqs.  (l)-(7)  are: 

At  the  ‘in’  terminal: 


j 


V 


oc 


MRi  +  Rl) 


and 


(8) 


where  Voc  is  the  Seebeck  voltage,  the  total  built-in  open  circuit  volt¬ 
age  in  the  iTED  at  no-load  condition.  In  Eq.  (8),  £,  A  and  RL  represent 
the  direction  normal  to  the  surface,  the  cross-sectional  area  and  the 
external  load  resistance,  respectively. 

The  Voc  is  calculated  as  the  summation  of  the  Seebeck  potentials 
at  the  interfaces  of  the  semiconductor  and  conductor  materials  and 
it  is  written  as: 


w  = 


2N  , 

Y.Y.X  l 

j-n,p  i-1  n>'  jAH 


2N,  1  f  .  .dT.. 
—  I  \<Xji\-^dA 


(9) 


The  Seebeck  potentials  are  evaluated  with  reference  to  the  hot 
surface  temperatures  at  the  interface  of  the  semiconductor  and  in¬ 
ter-connector  materials.  In  Eq.  (9),  A  is  the  interface  surface  area  at 
the  junction  between  the  thermoelectric  material  and  the 
connector. 

Further,  in  Eq.  (8),  Rt  is  the  total  internal  resistance  which  is  the 
summation  of  the  electrical  resistivities  of  the  n-  and  p- type  slices, 
the  connectors  and  the  inter-connector  materials  and  it  is  evalu¬ 
ated  as: 


2  N 


*=  E  E 

j=n,p,c,ic  i=  1 


Hji 

4 


i  [  Padv 

vji  Jv» 


(10) 


where  H,  A  and  V  are  the  height,  cross-sectional  area,  and  volume  of 
the  electrical  material,  respectively.  In  Eqs.  (9)  and  (10),  j  is  the 
material  type  (n,  p,  connector  and  inter-connector),  i  is  the  number 
of  vertical  legs  and  N  is  the  number  of  modules,  i.e.  each  module  has 
two  vertical  legs  as  shown  in  Fig.  la. 

At  the  ‘out’  terminal: 


V  =  0  and  ^-  =  0  (11) 

At  the  top  surface  of  the  upper  connector  and  the  bottom  sur¬ 
faces  of  the  lower  connectors: 


T  =  TC  and  g«0  (12) 

At  all  other  surfaces  exposed  to  surroundings: 

^  =  0  and  ^  =  O(iTED);  ^  =  0  main  flow  channel  (13) 
At  the  fluid  flow  inlet: 


x  =  0  :  u  =  U,  v  =  w  =  0 


and 


T  =  Tin 


(14) 


At  the  fluid  flow  outlet: 


x  =  5 W  +  N(2W  +  L)  +  20W  :  ^7  =  ^=v  =  w  =  0  (15) 

At  the  interface  between  the  semiconductor  and  connector  or  in¬ 
ter-connector  materials,  the  continuity  of  temperature,  current 
density  and  the  heat  flux  conditions  are  imposed  and  are  written  as: 


Tc,ic  —  Tn,p ,  Jc,ic  —  Jn,p 


dTc,ic 

d( 


kn.p  dTn  p 
Kjc~df 


(16) 


The  power  output  P0  from  the  iTED  for  a  given  load  resistance  RL 
and  the  heat  input  Qh  into  the  iTED  are  evaluated  as: 


(17) 


where  As  is  the  surface  area  of  inter-connector  walls  exposed  to  hot 
fluid. 

The  thermoelectric  conversion  efficiency  of  an  iTED  is  calcu¬ 
lated  as 


(18) 


3.  Numerical  solution  procedure,  grid  independence  study,  and 
validation 

The  numerical  simulations  are  performed  using  the  finite  vol¬ 
ume  formulation  of  Eqs.  (l)-(6)  and  the  constitutive  relation  (7) 
along  with  associated  boundary  conditions  Eq.  (8)  and  Eqs.  (11)- 
(15)  in  the  FLUENT-UDS  (User  Defined  Scalar)  environment.  The 
motivation  for  implementing  the  present  numerical  code  in  FLU¬ 
ENT  is  because  this  package  has  already  well-defined  fluid  and  heat 
transfer  models,  hence  it  allows  the  detailed  investigation  of 
fluid-thermo-electrical  coupled  phenomenon  in  iTEDs.  The  pres¬ 
sure-velocity  coupling  is  handled  using  SIMPLE  algorithm  [26]. 
The  spatial  descretisation  of  convective  and  diffusive  terms  is  done 
with  a  power  law  scheme  and  the  pressure  term  is  handled  with  a 
standard  scheme.  The  geometric  models  and  mesh  are  generated  in 
Gambit  2.4.  The  Seebeck  potential  and  Ohmic  electric  potential  dis¬ 
tributions  (Eq.  (7))  and  continuity  of  current  density  (Eq.  (2))  are 
calculated  using  UDS  fields.  The  Ohmic  heating,  Peltier  and  Thom¬ 
son  effects  are  modeled  as  source  terms  in  the  energy  equation  (Eq. 
(6)).  The  electric  current  is  evaluated  based  on  the  Seebeck  voltage 
(Eq.  (9))  produced  at  a  given  load  resistance  RL,  as  given  in  Eq.  (10). 
A  further  detailed  description  on  numerical  implementation  of  TED 
governing  equations  in  the  FLUENT-UDS  has  been  given  in  the  arti¬ 
cles  [17,19]. 

The  convergence  criteria  (the  difference  between  consecutive 
iterations  over  a  domain)  for  mass  and  momentum,  energy,  current 
density,  Seebeck  and  Ohmic  electric  potentials  are  set  as  10-5, 
10-15, 10-10, 10-10,  and  10-10,  respectively.  After  performing  a  grid 
independence  study,  suitable  grid  sizes  were  chosen  for  further 
numerical  simulations  on  iTEDs.  Grid  sizes  with  1,574,424, 
2,303,208,  2,906,712,  3,760,776  and  4,489,560  cells  are  used  for 
one-to-five  modules  iTEDs,  respectively,  to  generate  results.  For 


Table  1 

Grid  independence  study  for  single-stage  integrated  thermoelectric  device  (grid  size  in  bold  face  is  chosen  for  further  simulations)  [Re  =  500,  Tin  =  550  K,  Tc  =  300  K,  d  =  5  mm  and 
Rl  =  6.5  x  10“3  QJ. 


Cells 

PoW 

%  of  error  in  P0 

Qh  w 

%  of  error  in 

YJ  % 

%  of  error  in  rj 

406,572 

0.0929 

2.8211 

3.2932 

713,880 

0.0949 

2.1072 

2.7833 

1.3577 

3.4098 

3.4185 

1,574,424 

0.0959 

1.0828 

2.7652 

0.6556 

3.4697 

1.7271 

2,392,952 

0.0955 

0.4244 

2.7594 

0.2074 

3.4622 

0.2166 
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Table  2 

Dimensions  of  an  integrated  thermoelectric  device. 


H  (mm) 

D  (mm) 

W  (mm) 

t  (mm) 

L  (mm) 

d  (mm) 

20 

5 

5 

1.5 

10 

0.25,  1,  3,  5,  7 

brevity,  the  grid  independence  test  conducted  on  a  single-stage 
iTED  (as  shown  in  Fig.  1  a)  is  depicted  in  Table  1 .  Further,  the  orthog¬ 
onal,  nonuniform  grid  used  for  a  single-stage  iTED  along  the 
sectional  plane  (as  shown  in  side  view  of  Fig.  la)  is  also  shown  in 
Fig.  lb.  The  implemented  thermoelectric  model  has  been  validated 
with  published  results  and  are  given  in  our  previous  works  [19,23]. 


4.  Results  and  discussion 

Using  the  implemented  fluid-thermo-electric  coupled  field 
numerical  code  in  the  FLUENT  UDS  environment  as  presented  in 
the  previous  section,  the  thermoelectric  performance  of  an  inte¬ 
grated  thermoelectric  device  (iTED)  applied  to  waste-heat  recovery 
has  been  investigated.  In  this  study,  the  effects  of  various  parame¬ 
ters:  hot  fluid  inlet  temperatures  350  ^  Tin{I<)  <  550,  semiconduc¬ 
tor  (n-  and  p- type)  element  sizes  0.25  <  d(mm)  ^  7  and  the 
number  of  modules  1  ^  N  ^  5  on  the  performance  of  an  iTED 
expressed  in  terms  of  power  output  P0,  heat  input  Q^,  thermal  con¬ 
version  efficiency  p,  pressure  drop  A P/L,  produced  electric  current  I 
and  the  Ohmic  and  Seebeck  voltages  V  have  been  investigated  for 
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Fig.  2.  Effects  of  Re  on  (a)  x-velocity,  (b)  temperature  contours,  (c)  Ohmic  and  (d)  Seebeck  voltage  distributions  [at  Tin  =  450  K,  Tc  =  300  K,  d  =  5  mm  and  RL  =  /?,■]. 
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Fig.  3.  Effects  of  Re  on  (a)  heat  input  and  load  resistance,  (b)  power  output  and  thermal  efficiency  and  (c)  electric  current  and  voltages. 
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Fig.  4.  Influence  of  hot  fluid  inlet  temperature  Tin  on  (a)  heat  input  and  load  resistance,  (b)  power  output  and  thermal  efficiency  and  (c)  electric  current  and  voltages. 


various  fluid  flow  configurations  50  <  Re  <  500.  The  dimensions  of 
an  iTED  used  in  the  numerical  simulations  (as  shown  in  Fig.  la)  are 
given  in  Table  2. 


The  thermoelectric  materials  are  treated  as  bulk  bismuth-tellu- 
ride-antimony-selenium:  n-type  (75%Bi2Te325%Bi2Se3)  andp-type 
(25%Bi2Te375%Sb2Te3  (1.75%  excess  Se)).  Copper  is  used  for  the 
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conductor  and  air  is  the  hot  fluid.  Further,  the  thermo-electrical 
properties  of  the  semiconductor  (Seebeck  coefficient  a,  electrical 
resistivity  p  and  thermal  conductivity  k)  and  conductor  (p  and  k) 
materials  are  temperature-dependent  (as  reported  in  article  [23], 
in  Table  1)  while  the  thermo-physical  properties  of  air  are  kept 
constant.  For  the  numerical  simulations,  the  maximum  power 
output  P0  from  the  iTED  is  achieved  when  the  load  resistance  RL 
equals  the  total  internal  resistance  Rt  (Eq.  (10)),  and  thus  all  cases 
use  this  relation. 

For  a  single-stage  iTED,  the  contours  of  x-velocity  and  isotherms 
at  section  z  =  2.5  mm  and  Ohmic  and  Seebeck  voltage  distributions 
at  z  =  0.5  mm  are  shown  in  Fig.  2  for  Re  values  of  1 00  and  500.  Here, 
Tin  =  450  K  and  d  =  5  mm  are  kept  constant.  For  clarity,  as  displayed 
in  Fig.  2a  and  b,  an  upstream  length  W  and  a  downstream  length 
5  W  are  chosen  for  depicting  x-velocity  and  temperature  distribu¬ 
tions.  From  Fig.  2b,  it  is  seen  that  an  increase  in  temperature  gra¬ 
dients  in  the  fluid  region  at  the  inter-connector  channel  walls 
when  Re  is  increased  from  100  to  500.  This  increase  in  temperature 
gradients  helps  in  enhancement  of  convective  heat  transfer  and 
thus  eventually  results  in  a  reduction  of  the  total  thermal  resis¬ 
tance  of  the  iTED.  Furthermore,  the  recirculation  zones  formed  in 
the  up  and  down  streams  of  iTED  and  in  between  the  legs  have 
minimal  effect  on  device  performance.  However,  the  strength  of 


the  recirculation  increases  with  an  increase  in  Re  values,  as  seen 

in  Fig.  2a. 

The  increase  in  temperature  differential  Th-Tc  (between  the 
inter-connector  channel  wall  and  the  cold  end  surface)  via  increas¬ 
ing  Re  generates  higher  Seebeck  voltages  in  the  thermoelectric 
elements,  as  shown  in  Fig.  2d.  In  each  element,  the  Seebeck  voltage 
drop  is  calculated  with  respect  to  the  connector  (copper)  material 
value;  however,  the  Seebeck  coefficient  of  copper  is  negligible 
compared  to  the  semiconductor  value.  From  Fig.  2d,  at  given  a 
Re,  it  is  noticed  that  due  to  presence  of  a  larger  Th-Tc  value  in  the 
left  leg,  the  built-in  Seebeck  voltage  is  higher  compared  to  the  right 
leg  value.  On  the  other  hand,  as  shown  in  Fig.  2c,  as  Re  increased 
from  100  to  500,  the  Ohmic  voltage  drop  increased  from  0.017  to 
0.031  V.  Further,  the  Ohmic  voltage  drop  in  the  copper  material 
is  negligible  due  to  its  lower  electrical  resistivity  value  (approxi¬ 
mately  three  orders  of  magnitude  smaller  than  semiconductor 
material  value). 

The  effects  of  several  flow  configurations  Re  on  iTED  perfor¬ 
mance  parameters  Qh,  Po  and  r\,  and  electrical  characteristics  I 
and  V  are  shown  in  Fig.  3.  The  change  in  total  internal  resistance 
Ri  and  the  pressure  drop  A P/L  with  Re  are  also  depicted  in  Fig.  3a 
and  b,  respectively.  Fig.  3a  represents  the  heat  transferred,  as  given 
by  Qhi  and  Qh2,  to  the  left  and  right  legs  of  a  single-stage  iTED 


(b) 


X- Velocity  (C)  Isotherms 


Fig.  5.  The  variation  ofx-velocity  and  temperature  contours  for  different  thermoelectric  element  heights  d  (mm)  (a)  0.25  (b)  3  and  (c)  7  of  a  single-stage  thermoelectric  device 
[at  Re  =  500,  Tin  =  550  K,  Tc  =  300  K  and  RL  =  /?,-]. 
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module,  respectively.  For  a  given  Tin  =  450  K,  an  increase  in  Re 
shows  an  enhancement  in  Qh  and  resulted  in  higher  Th-Tc  values. 
The  increase  in  Th-Tc  is  due  to  the  reduction  of  developing  ther¬ 
mal  boundary  layer  thickness  over  inter-connector  channel  walls 
as  Re  increases  (as  shown  in  Fig.  2b).  As  discussed  earlier,  the 
higher  Th-Tc  values  result  in  larger  built-in  Seebeck  voltages 
which  in  turn  results  in  a  greater  P0  (as  seen  in  Fig.  3b).  At 
Re  =  500,  P0  and  r\  increase  five-  and  twofold,  respectively,  as 
compared  to  the  values  at  Pe  =  100.  The  produced  current  I 
and  voltages  V  have  shown  similar  trends  with  Re.  From 
Fig.  3c,  it  is  noticed  that  at  a  given  Re  value,  the  Ohmic  voltage 
drop  is  higher  than  the  Seebeck  voltage  drop  value.  Furthermore, 


the  flow  rate  has  a  minimal  effect  on  Rt  predictions  and  merely  a 
3.3%  variation  is  observed  in  Rt  for  the  flow  rate  range  of 
50  ^  Re  <  500.  This  change  merely  is  due  to  the  variation  in  elec¬ 
trical  resistivity  value  with  temperature  (Eq.  (10)).  As  seen  in 
Fig.  3b,  the  A P/L  increases  as  Re  is  increases  and  shows  a  nonlin¬ 
ear  behavior.  The  reason  for  this  behavior  is  due  to  the  cumula¬ 
tive  effect  of  developing  flow  in  the  integrated  flow  channels  and 
the  formation  of  the  recirculation  zones  in  the  main  flow  chan¬ 
nel  (as  shown  in  Fig.  2a). 

For  a  given  RL  =  Rt  condition,  the  variation  in  both  thermal  and 
electrical  characteristics  of  a  single-stage  iTED  module  with  hot 
fluid  inlet  temperatures  Tin  are  shown  in  Fig.  4a-c.  An  increase  in 
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Fig.  6.  The  variation  of  Ohmic  and  Seebeck  voltage  distributions  for  different  thermoelectric  element  heights  d  (mm)  (a)  0.25  (b)  3  and  (c)  7  of  a  single-stage  thermoelectric 
device  [at  Re  =  500,  Tin  =  550  K,  Tc  =  300  K  and  RL  =  Rt ]. 
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Tin  results  in  an  enhancement  in  Qh,  Po,  *7,  I  and  V  values.  Further,  a 
linear  behavior  in  V-I( Eqs.  (8)),  Qh  (Eq.  (17))  and  r\  (Eq.  (18))  values 
and  a  nonlinear  trend  in  P0  (Eq.  (17))  and  ft,  (Eq.  (10))  with  Tin  are 
observed.  The  reason  for  this  nonlinear  trend  in  ft,  is  due 


to  the  behavior  of  temperature-dependent  electrical  resistivity. 
Moreover,  the  increment  in  ft,  with  an  increase  in  Tin  is  marginal 
and  exhibited  a  9.2%  variation  within  the  temperature  range 
350  ^  !,„(!<)  <  550  (as  seen  in  Fig.4a). 


d,  mm 
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Fig.  7.  The  response  of  (a)  heat  input  and  load  resistance,  (b)  power  output  and  thermal  efficiency  and  (c)  electric  current  and  voltages  for  different  thermoelectric  element 
heights  d. 
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Fig.  8.  The  influence  of  number  of  modules  N  on  (a)  heat  input  and  load  resistance,  (b)  power  output  and  thermal  efficiency  and  (c)  electric  current  and  voltages. 
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Figs.  5  and  6  show  the  contours  of  x-velocity  and  temperatures 
at  section  z  =  2.5  mm  and  the  Ohmic  and  Seebeck  potential  distri¬ 
butions  at  z  =  0.5  mm,  respectively  for  three  sizes  of  semiconductor 
elements  d.  Here,  both  the  n-  and  p- type  semiconductor  elements 
are  kept  equal  in  size  and  also  the  height  of  the  module  is  main¬ 
tained  constant.  Hence,  the  modification  in  d  effects  inter-connec¬ 
tor  channel  dimension  (H  -  2d)  xWxD  as  shown  in  Fig.  la.  For 
fixed  Tin  =  550  K  and  Re  =  500,  where  Re  is  defined  based  on 
hydraulic  diameter  of  the  main  flow  channel  ( H  -  2d)  x  D,  as  d  in¬ 
creases,  both  averaged  velocity  of  fluid  in  the  inter-connector 
channel  and  total  thermal  resistance  of  legs  increase  and  hence, 
an  increase  in  the  inter-connector  channel  wall  temperature  is  no¬ 
ticed.  Further,  the  gain  in  Th-Tc  with  an  increase  in  d  helps  to  en¬ 
hance  the  total  built-in  Seebeck  potential  in  the  thermoelectric 
materials  and  also  to  achieve  higher  Ohmic  voltage  drop  values 
(as  shown  in  Fig.  6).  It  is  noticed  that  an  increase  in  d  changes 
the  developing  flow  behavior  in  the  inter-connector  channel,  and 
thus  the  temperature  distribution  in  the  fluid  and  solid  regions  ap¬ 
pears  accordingly.  The  strength  of  recirculation  in  between  the  legs 
and  the  reattachment  length  in  the  downstream  decrease  with  an 
increase  in  d  value. 

The  response  of  iTED’s  thermoelectric-hydraulic  characteristics 
for  various  thermoelectric  material  sizes  d  are  presented  in  Fig.  7.  As 
seen  in  Fig.  7a,  the  increase  in  d  results  in  a  linearly  proportional 
gain  in  ft*  values.  Further,  this  Rt  has  a  phenomenal  effect  on  iTED’s 
performance.  When  Re  and  Tin  are  fixed,  the  Qh  decreases  with  an 
increase  in  d  which  is  due  to  the  increase  in  the  thermal  resistance 
of  the  legs.  However,  at  a  given  d,  the  left  leg  extracts  more  heat 
compared  to  right  leg  of  the  iTED.  From  Fig.  7b,  it  is  observed  that 
there  exists  a  maximum  P0  at  an  optimum  d  value  for  a  given  set 
of  operating  parameters.  This  behavior  is  due  to  the  mutual  inter¬ 
play  of  the  total  electrical  resistance  and  the  thermal  resistance  val¬ 
ues  of  an  iTED.  Here,  the  total  thermal  resistance  includes  the 
summation  of  resistances  due  to  semiconductor  elements,  connec¬ 
tors,  convective  heat  transfer  at  the  inter-connector  channel  walls 
and  the  constant  temperature  at  the  cold  surfaces.  On  the  other 
hand,  The  r\  of  the  iTED  increases  exponentially  with  an  increase 
in  d.  For  a  fixed  Re,  the  effect  of  d  on  AP/L  is  marginal  and  it  is 
due  to  the  increment  in  inlet  velocity  and  the  reduction  in  hydraulic 
diameter  with  an  increase  in  d  value.  The  Th-Tc  value  improves  with 
an  increase  in  d  and  thus  reflects  in  the  total  Seebeck  voltage  poten¬ 
tial  predictions  (Fig.  7c).  As  shown  in  Fig.  7c,  the  produced  current  / 
reduces  and  the  Ohmic  voltage  potential  increases  with  a  rise  in  d 
values  and  these  trends  are  due  to  the  increment  in  ft,  values. 

The  performance  of  a  multi-stage  iTED  for  a  various  number  of 
modules  (1  ^  N  ^  5)  connected  electrically  in  series  and  thermally 
in  parallel  is  shown  in  Fig.  8.  For  all  other  parameters  kept  invari¬ 
ant  (Re  =  500,  Tin  =  550  K,  Tc  =  300  K  and  d  =  5  mm),  the  ft,  of  the 
iTED  increases  linearly  with  the  addition  of  N  and  shows  a  linear 
behavior.  From  Fig.  8a,  it  is  noticed  that  the  change  in  heat  transfer 
rates  (Q/n-Q/is)  at  the  inter-connector  walls  of  each  leg  is  minimal 
with  the  addition  of  N.  However,  the  total  Qh  increases  with  the 
addition  of  N.  The  P0,  Ohmic  and  Seebeck  potentials  increase  signif¬ 
icantly  and  illustrate  more  than  two,  three  and  threefold  increases, 
respectively,  at  N  =  5  compared  to  N  =  1  values.  On  the  other  hand, 
both  the  t]  and  the  /  decrease  by  nearly  33%,  at  N  =  5  when  com¬ 
pared  to  N  =  1  values.  Furthermore,  a  linear  trend  in  AP/L  with  an 
addition  of  modules  N  is  observed. 


5.  Conclusions 

The  performance  of  an  integrated  thermoelectric  device  (iTED) 
applied  to  waste-heat  recovery  has  been  investigated  using  three- 
dimensional  fluid-thermo-electric  coupled  field  numerical  simula¬ 
tions.  The  influence  of  hot  fluid  temperatures  Tin,  thermoelectric 


elements  sizes  d  and  the  number  of  modules  N  on  iTED’s 
performance  have  been  studied  for  different  fluid  flow 
configurations  Re. 

An  increment  in  either  Re  or  Tin  shows  an  enhancement  in  the 
device  performance.  The  increase  in  Re  or  Tin  produces  larger 
temperature  differentials  (Th-Tc)  in  thermoelectric  materials  and 
thus  resulted  in  higher  Seebeck  voltages  which  in  turn  produced 
greater  power  output.  Further,  a  linear  behavior  in  produced 
electric  current  J,  voltages  V  and  efficiency  r\  values  and  a  non¬ 
linear  trend  in  power  output  P0  and  heat  input  Qh  with  Tin  are 
observed. 

An  increase  in  d  results  in  an  increase  in  both  thermo-electric 
resistances  and  Seebeck  voltages  and  a  decrease  in  Qh  values.  For 
a  given  set  of  conditions,  it  is  noticed  that  there  exists  an  optimum 
d  value  where  a  maximum  P0  occurs.  An  exponentially  increase 
trend  in  r\  and  voltages  and  a  linear  decrease  in  I  with  an  increase 
in  d  are  observed.  Further,  the  effect  of  d  on  the  pressure  drop  AP/L 
is  shown  marginal. 

As  the  number  of  modules  N  increases,  ft„  Qh,  Po  and  the  Ohmic 
and  Seebeck  voltages  increase  substantially;  conversely,  r\  and  /  de¬ 
crease  as  a  response  to  the  increasing  aforementioned  values. 

It  is  recommended  to  use  fluid  flow,  heat  transport  and  electri¬ 
cal  coupling  field  simulations  to  capture  the  thermoelectric  perfor¬ 
mance  of  state-of-the-art  integrated  thermoelectric  devices. 
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